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Abstract. This article is concerned with an oceanographic model de- 
scribing the asymptotic behaviour of a rapidly rotating and incompress- 
ible fluid with an inhomogeneous rotation vector; the motion takes place 
in a thin layer. We first exhibit a stationary solution of the system 
which consists of an interior part and a boundary layer part. The spa- 
tial variations of the rotation vector generate strong singularities within 
the boundary layer, which have repercussions on the interior part of 
the solution. The second part of the article is devoted to the analysis 
of two-dimensional and three-dimensional waves. It is shown that the 
thin layer effect modifies the propagation of three-dimensional Poincare 
waves by creating small scales. Using tools of semi- classical analysis, 
we prove that the energy propagates at speeds of order one, i.e. much 
slower than in traditional rotating fluid models. 

Resume. On etudie ici le comportement asymptotique d'un fluide incom- 
pressible tournant a grande vitesse dans une couche mince, avec un vec- 
teur rotation inhomogene ; ce type de modele apparait en oceanographie. 
On commence par exhiber une solution stationnaire du systeme, obtenue 
comme la somme d'un terme interieur et d'un terme de couche limitc. 
Les variations spatiales du vecteur rotation generent de fortes singula- 
rites dans la couche limite, qui se repercutent dans la partie interieure 
de la solution. Dans un second temps, on caracterise le comportement 
des ondes bi- et tri-dimensionnelles. L'effet de couche mince modifie la 
propagation des ondes de Poincare (3D) en favorisant l'apparition de 
petites echelles. Grace a une analyse de type semi-classique, on montre 
que la vitesse de propagation de l'energie est d'ordre un, soit beaucoup 
plus faible que dans les modeles classiques de fluides tournants. 

1. Introduction 

The goal of this article is to study the behaviour of a rotating, incom- 
pressible and homogeneous fluid, whose rotation vector depends on the 
(horizontal) space variable. We also assume that the motion of the fluid 
takes place in a thin layer. These two features are inspired from models of 
oceanic circulation, which are the main physical motivation for our study. 
We will explain more thoroughly the physical assumptions and scalings lead- 
ing to our model in paragraph 11.11 
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The mathematical framework of our analysis is the following: consider 
the equation 



{x h ,z) 6U(,x (0,1), 

where the horizontal domain ujh is either T 2 or T x R. Equation (|1.1[) is 
endowed with Navier conditions at the bottom of the domain 



and we assume that there is a shear stress at the surface of the fluid, de- 
scribed by the boundary condition 



Above, e, rj, Uh, u z , 7 are positive parameters, whose relative size will be pre- 
cised later on. Let us merely announce that e, rj, Uh, v z are meant to be small, 
whereas 7 will be taken large. We emphasize that equation (|l.ip . supple- 
mented with (jl.3p - (ll.2p . is already in rescaled form. Hence all quantities are 
dimensionless. We refer to the next subsection for a derivation of this equa- 
tion, and for a definition of the various parameters in terms of the physical 
quantities involved in the model. 

Notice that the rotation is of order e , with t <1; hence we focus on 
the limit of high rotation. As we will see in paragraph ll.lt the parameter 
rj is the aspect ratio of the domain: assuming that rj <C 1 means that the 
characteristic horizontal length scale is much larger than the vertical one. 
In other words, the motion is set in a thin layer. 

In this article, we are primarily interested in two topics: the computation 
of stationary solutions of our model, and the analysis of the local stability 
of these stationary solutions in the case coh = T x R. In particular, we will 
not address the full Cauchy problem here. Indeed, it can be proved that in 
the scaling which is the most relevant for our study, the energy estimates 
for the system (jl.ip - (jl.2j) - (jl.3p explode in finite time. In a similar way, the 
stationary solution that we build has a size which becomes arbitrarily large 
as e,rj vanish. Hence the problem (|l.ip - (jl.2p - (|1.3p is highly singular. 

To our knowledge, the asymptotic analysis of the system (|l.ip has not 
been addressed before: in the papers [10] by I. Gallagher and the second 
author, and then [8] by A. Dutrifoy, A. Majda and S. Schochet, the authors 
study the asymptotic behaviour of a shallow water system within a /3-plane 
model (i.e. in the case b(xh) = f^x^). This shallow water system can be 
obtained by considering the limit rj — ► in (II. ip (see [H]). Thus the studies 
of |10} [8] are concerned with the successive limits r\ — > 0, e — > 0. In [7], B. 
Desjardins and E. Grenier take into account the thin layer effect within the 
original Navier-Stokes system, but they assume that b{xh) = 1 + ex2] hence 



(1.1) 




(1.2) 



9zUh\z=0 = °> u 3\z=0 = 



(1.3) 



d z u h \ z= i(t,x h ) = ja(xh), 
u 3 \z=i = 0- 
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the penalization is constant at first order. Our goal is to study a crossed 
limit (e, rj) — > (0, 0), with a rotation vector which has variations at the main 
order. 

Let us now make precise the main novelties of our work: first, the con- 
struction of stationary solutions involves the definition of boundary layer 
terms with a varying Coriolis factor b. Since the size of the bound- 
ary layer is directly related to the amplitude of b, singularities appear at 
the vanishing points of b. These singularities in the boundary layer have 
repercussions on the interior part of the stationary solution, and make the 
construction much more involved than in the constant case. On the other 
hand, studying the stability of stationary solution when w/, = TxR amounts 
to describing the waves in the /3-plane model with a thin layer effect. 
We exhibit new types of behaviour for the Poincare waves, for which we 
prove that dispersion takes place on a time scale much larger than usual: 
for instance, in Chapter 2 of [11], the group velocity associated with Poincare 
waves (i.e. the speed at which energy propagates) is of order e _1 , while the 
group velocity in the present setting is of order one. The proof of this fact 
uses tools of semi-classical analysis, in the spirit of the recent papers by C. 
Cheverry, I. Gallagher, T. Paul and the second author (see [5j[6]). Notice 
also that the presence of dispersion in an oceanographic model is itself un- 
usual: indeed, most models are set in a compact domain (see [3]), where 
no dispersion can occur. Moreover, the most commonly used whole-space 
model is the shallow water system within the /3-plane model (see |10t [8]), 
for which waves are trapped into a waveguide, and thus no dispersion occurs 
either. 

In the next paragraphs, we explain which physical assumptions led to the 
system (11. ip . We then present our main results. Eventually, let us point 
out that the structure of the stationary solution which will be built in this 
article enforces particular shapes for the isothermal surfaces inside the fluid 
(the so-called "thermocline" ) . We present a few results in this regard in 
paragraph 11.31 

1.1. Physical derivation of equation (II. ip . 

Let us now explain in which regime oceanic currents can be modeled by 
equation (II. ip . In this subsection, we denote by u the velocity of oceanic 
currents in dimensional variables. The dimensionless variables, i.e. the ones 
in which equation (11.11) is written, will be denoted with a prime. 

• As a starting point, we recall that the ocean can be considered as an 
incompressible fluid with variable density p. In order to simplify the analysis, 
we neglect the variations of density, which are of order 10~ 3 in the ocean. 
Consequently, the velocity u satisfies the Navier-Stokes equations, with a 
Coriolis term accounting for the rotation of the Earth 



(1.4) 



po [d t u + (u ■ X7)u] + Vp = T + p u A Q , 

V -u = 0, 
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where T denotes the frictional force acting on the fluid, O is the (vertical 
component of the) Earth rotation vector, p is the pressure defined as the 
Lagrange multiplier associated with the incompressibility constraint, and po 
is the (constant) value of the density. 

Since we have chosen to work on large horizontal scales (see below), equa- 
tion (11. 4p should be written in spherical coordinates. However, computations 
involving spherical coordinates are much lengthier, and do not change sub- 
stantially the physical phenomena we wish to highlight, at least at a formal 
level (see [23 1). Thus in the rest of the article, we neglect the curvature 
of the Earth (but we keep a varying Coriolis factor nonetheless). Note 
also that we neglect the influence of the horizontal component of the Earth 
rotation vector, which is classical in an oceanographic framework (see |11|). 

The observed persistence over several days of large-scale waves in the 
oceans shows that frictional forces T are weak, almost everywhere, when 
compared with the Coriolis acceleration and the pressure gradient, but large 
when compared with the kinematic viscous dissipation of water. One com- 
mon but not very precise notion is that small-scale motions, which appear 
sporadic or on longer time scales, act to smooth and mix properties on the 
larger scales by processes analogous to molecular, diffusive transports. For 
the present purposes it is only necessary to note that one way to estimate 
the dissipative influence of smaller-scale motions is to retain the same rep- 
resentation of the frictional force 

T = A h A h u + A z d zz u 

where A z and Ah are respectively the vertical and horizontal turbulent vis- 
cosities, of much larger magnitude than the molecular value, supposedly 
because of the greater efficiency of momentum transport by macroscopic 
chunks of fluid. Notice that A z ^ Ah is therefore natural in a geophysical 
framework (see |23j). Moreover, models of oceanic circulation usually as- 
sume that the vertical viscosity A z is not constant (see [21 [22] ) ; we choose to 
retain only the mean boundary value of the vertical viscosity A z , since one 
of the motivations for our work was to compute the boundary layer terms 
in a context where Q is not constant. 

• Let us now describe the boundary conditions associated with (jl.4p : 
typically, Dirichlet boundary conditions are enforced at the bottom of the 
ocean and on the lateral boundaries of the horizontal domain ujh (the coasts), 
i.e. 

^ u ]z=hB{xh) =0 (bottom), 

u \xedw h = (coasts). 

In equation we have neglected the effects of the lateral boundary 
conditions by considering the case when Uh is either T x R or T 2 . By 
doing so, we have deliberately prohibited the apparition of strong western 
boundary currents, which play a crucial role in the oceanic circulation (e.g. 
the Gulf Stream, the Kuroshio current). These horizontal boundary layers 
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are believed to be responsible for the vertical structure of the ocean, and for 
the creation of large eddies. In the linear case, the mathematical treatment 
of these layers, called Munk layers, is performed by B. Desjardins and E. 
Grenier in [7]. Their study could probably be mimicked in the present paper 
without strong modifications; however, we have chosen to leave this issue 
aside in order to focus on the other features of the model. Note that in the 
nonlinear case, the analysis of lateral boundary layers is completely open 
from a mathematical point of view. 

In a similar fashion, for the sake of simplicity, we did not take into account 
the topography of the bottom in (|1.2p (i.e. we have taken hs = 0), 
and we took Navier instead of Dirichlet boundary conditions, meaning that 
oceanic currents achieve perfect slip on the bottom. This choice simplifies 
the mathematical analysis, since it avoids the apparition of Ekman boundary 
layers on the lower boundary. The treatment of Ekman boundary layers in 
the case of a Dirichlet boundary condition with hs = is in fact completely 
similar to the one of Ekman boundary layers due to the wind at the surface 
of the fluid, which is performed in section [2 Hence changing Dirichlet into 
Navier boundary conditions is not a strong mathematical restriction. The 
case of Ekman boundary layers with a non-zero Hb has been addressed by B. 
Desjardins and E. Grenier [7], N. Masmoudi [20], and D. Gerard- Varet [13] 
in the case of a constant b, when \ib is of the order of the Ekman boundary 
layer (see below). In the present case, if the same assumption is satisfied, it 
can be checked that the case of a non-constant hs can be treated with the 
same arguments as the ones in section [2j 

We assume that the upper surface, which we denote by T s , has an equation 
of the type z = hs(t,Xh)- As boundary conditions on T s , we enforce (see 

EH!) 

S • nr s = a w , 

-Q^-0<z<h s (t,x) + ^ v x{i-0<z<h s (t,x h ) u ) = U 

where X is the total stress tensor of the fluid, and cr w is a given stress tensor 
describing the wind on the surface of the ocean. In general, T s is a free 
surface, and a moving interface between air and water, which has its own 
self consistent motion. In (11.31). we have assumed that 



hs(t,x h ) = D, 

where D is the typical depth of the ocean. Hence (|1.3j) is a rigid lid approx- 
imation, which is a drastic, but standard simplification. The justification of 
(jl.3p starting from a free surface is mainly open from a mathematical point 
of view; we refer to [l] for the derivation of Navier-type wall laws for the 
Laplace equation, under general assumptions on the interface, and to [16] for 
some elements of justification in the case of the great lake equations. Nev- 
ertheless, from a physical point of view, the simplification does not appear 
so dramatic, since in any case the free surface is so turbulent with waves 
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and foam, that only modelization is tractable and meaningful. Condition 
(jl.3p is a simple modelization which already catches most of the physical 
phenomena (see |23j). 

• Let us now evaluate the order of magnitude of the different parameters 
occurring in (II. 4ft . and write the equations in a dimensionless form. We set 

Uh = Uu' h , u 3 = Wu' 3 , 

Xh = Hx' h , z = Dz', 

where U (resp. W) is the typical value of the horizontal (resp. vertical) 
velocity, H is the horizontal length scale, and D the depth of the ocean. In 
order that u'{x') remains divergence-free, we choose 

UD 

w = —. 

A typical value of the horizontal velocity for the mesoscale eddies that have 
been observed in western Atlantic (see for instance [23]) is U ~ 1 cm • s _1 . 
Moreover, the typical horizontal and vertical scales which we are interested 
in are 

H ~ 10 4 km, and D ~ 4 km. 

Notice that we work on an almost planetary scale, which justifies the use of 
a varying rotation vector. Concerning the rotation, we write Vt = VLqsvo.(9), 
where 9 is the latitude, and Qq = 27r/day ~ 7 • 10~ 5 s . Eventually, we 
consider the motion on a typical time scale T, with T of the order of a few 
months (T ~ 10 7 s). With these values, we get 

1 



10" 



and hence e<l (notice that the parameter e is dimensionless). Thus the 
asymptotic of fast rotation (small Rossby number) is valid. 

Thus the dimensionless system (see for instance |23| 115]) becomes 
(1.7) 

d t u + ^jj-u ■ Vu + ^6(x/0e 3 A u' + (J^j^j/) ~ v h^hU - v z d zz u = 0, 

V • u' = 0, 

where rj := D / H ~ 4-10 is the aspect ratio, and the vertical and horizontal 
viscosities are defined by 

TA Z A h T 



Typical values for the turbulent viscosities are (see [E]) A z /po ~ 10~ 4 — 
10~ 3 m 2 • s —1 , and A^/po ~ 10 4 — 10 5 m 2 • s~ , which yields in the present 
case v z ~ 10 -3 and u h ~ 10~ 10 - 1CT 9 . 

The boundary conditions are (|1.3p . (|1.2[) . with 

_ \a w \D 
7 : ~ AM ' 
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Notice that with the time scale chosen above, the convective term is of order 
1CP 2 -C 1; hence we neglect it in the rest of the study. Note however that the 
effect of this term is expected to be large if the waves associated with (|1.7j) 
are resonant, and small if they are dispersive. Thus the rigorous treatment 
of the convective term requires a mathematical analysis which goes beyond 
the scope of this article, and which we deliberately leave aside from now on. 

In the rest of the article, the relative size of the parameters will be chosen 
as follows: the most important feature of our analysis is that r) and e are 
chosen of the same order. In order to keep the number of different small 
parameters to a minimum, we also choose to take v z = e, and 7 = e -2 ; with 
this last choice, the interior part of the stationary solution built in the next 
sections will be of order one. Concerning the size of u^, our analysis allows 
us to consider horizontal viscosities Vh = o(e), which is compatible with the 
orders of magnitude given above. 

1.2. Main results. 

We present here two types of results: first, we build an approximate "sta- 
tionary" solution of the system (jl.ip . endowed with the boundary conditions 
(|1.3p - (jl.2|) . The problem studied is rather different from the Cauchy prob- 
lem, since no initial data is prescribed. The goal is merely to compute a 
solution of (jl.ip . and to investigate its asymptotic behaviour as e vanishes. 

Once the behaviour of the stationary solution is understood, we study its 
local stability; since equation (jl.ip is linear, this is equivalent to studying the 
Cauchy problem for equation (jl.ip . with homogeneous Navier conditions at 
z = and z = 1. We then exhibit Rossby waves, which are essentially two- 
dimensional, and Poincare waves, which are fluctuations around the three 
dimensional part of the initial data, and which take place on a much larger 
time scale. 

Let us now state our result about stationary solutions: since the vertical 
viscosity is small (we take v z = e <C 1), it disappears from the asymptotic 
system. As a consequence, solutions of the limit system cannot satisfy the 
boundary conditions. Thus boundary layer terms are introduced, which 
restore the correct boundary conditions. Hence the stationary solution built 
here is composed of an interior part and a boundary layer part. 

We state our result in the case u>h = T x R, and explain below the 
Theorem the main differences when LOh = T 2 . Throughout the paper, we 
set 

u := uj h x (0, 1). 

Theorem 1.1 (Stationary solutions of (jl.7p ). Let uj^ = T x R. 

Assume that = o(e) and that rj = v z = e, 7 = e~ 2 . 
Let a £ H 2 (uj h ) n W^°°(uj h ) such that 

^ 18 ^ W(x,y)\, \d x a(x,y)\ <Cy 2 V(x,y)eu h , 

\d y a(x,y)\ <C\y\ V(x,y) E u h 



8 



ANNE-LAURE DALIBARD AND LAURE SAINT-RAYMOND 



and such that the following compatibility condition is satisfied 
(1.9) / a 1 (x,y)dx = Vy. 



Assume that the Coriolis factor b satisfies the following assumptions: 

b(x, y) = b(y) V(x, y) E u h , with b G W to ' c °°(R), 

(1.10) 6(y) ^ for y ^ 0, and 3C > 0, |6(y)| > C for \y\ > 1, 

3c > 0, c" 1 < b'{y) < c Vy, 6(y) ~ [5y for y -» 0. 

T/ien i/iere exists stationary functions (u stat ,p stat ) g L 2 (uj) n H (oS), such 
that u stat satisfies ((131) . (Ol) and 



^(y)«°< ) x + V,,p^ - e^ttf* - ^Arf* = r£ + r 2 

- ed 22 < at - i/ftAfcuf" = r£ + r| 5 

r£ = o(l) i nL 2 (u,), r 2 = (^) *n L 2 ([0, 1], tf" 1 K)), 
r 3 x = ofe" 1 ) in L 2 H, r 2 = o^ 1 ^) m L 2 ([0, 1], tf" 1 (u h )). 
Moreover, u stat can be decomposed as 

u stat = U BL + u tot } 

where u BL is a term located in a boundary layer of size e, in the vicinity of 
the surface, and u mt is an interior term. The functions u BL and u tnt satisfy 
the following estimates 

IK n *IU 2 (w) < C || a \\m(u h ), 

ii B L 1 1 ii ii 

(1-11) 11% IIl 2 (u/) < —rWWwtuh), 

V e 

ll M f L ||L 2 H < C|klliri(w fc )- 

If oJh = T 2 , the result remains true under slightly different conditions on 
a and b. More precisely, we assume that a E F 2 (T 2 ) satisfies (|TS|I . (Ol) . 
and that 

d(suppa, (T, 1/2)) > 

In other words, a vanishes in a neighbourhood of (x, 1/2) for all x £ T (and 
by periodicity, in a neighbourhood of (x, —1/2) also). 

We assume furthermore that b(x, y) = b{y) with 
(1.12) 

b £ L°°(T) and b £ W 2 '°°(K) MK C T compact s.t. d(K, 1/2) > 0, 
b{y) / for y / 0, and 3C > 0, |6(y)| > C for |y| > 1/4, 

o(y) ~ /% for y -> 0, 

Vif C T compact s.t. d(if, 1/2) > 0, 3c A - > 0, c^ 1 < 6'(y) < c K My G c X - 
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In other words, we do not assume that b S Ty 2, °°(T): b may have a 
discontinuity at y = 1/2. But we require that a vanishes in a neighbourhood 
of that singularity, so that all terms of the type o~b, cr/b, a/b' are well-defined 
and T 2 -periodic. 

Remark 1.2. (i) The assumptions (I1.10p -( fl.l2p on the Coriolis factor b 
is satisfied in two particular cases: 

• b{y) = (3y, with ujh = T x R: this approximation is particularly 
relevant for the motion of equatorial currents, and is used in par- 
ticular in [10], [8]. 

• b(v) = sin(7ry/2), with u>h = T 2 : this is the case of a real ocean, 
whose study takes place on a planetary scale. Of course, in this 
case, the effect of the curvature of the Earth should be taken into 
account, which we have chosen not to do here (see the discussion 
in the previous paragraph). 

(ii) Notice that in the above Theorem, it is assumed that the surface stress 
vanishes near y = 0. Although this assumption stems from mathemat- 
ical considerations, it is in fact quite reasonable in an oceanographic 
context. Indeed, it is a well-known phenomena that there are no steady 
surface winds near the equator: as trade winds coming from the North 
and South meet, they are heated and produce upward winds. The area 
of calm in the vicinity of the equator is called the Doldrums. 

(iii) The compatibility condition (jl,9p means that there is no zonal average 
wind. This condition is of course not realistic from a physical point 
of view, but it is the price to pay for working with a domain with no 
boundary in x. If the horizontal domain oj^ is replaced by [0, 1] X T or 
[0,1] x R, this condition disappears; the (mathematical) counterpart 
lies in the construction of the horizontal boundary layer terms, the 
so-called Munk layers discussed in the previous paragraph. 

(iv) In general, the size of the boundary layer term u BL is much larger 
than that of the interior term. This means that the greatest part of the 
energy is concentrated in a boundary layer located in the vicinity of the 
surface. In the original variables, it can be checked that the boundary 
layer carries an energy of order pqU 2 H^ , while the energy contained 
in the interior of the domain is of order pqU 2 H 2 D. 

This is in fact a consequence on the requirements on u BL , u mt , 
and not an artefact of our model. Indeed, assume that the functions 
u BL , u mt are such that 

1 1 ZTVt II || ITht || || ITht 1 1 

\\U3\z=l\\L2(u, h ) ~ \\U h \ z=1 \\L2(u, h ) ~ \\ u h IIl 2 (u/), 

int __ _„.BL 
u 3|z=l ~~ u 3|.z=l' 

and assume that u BL ,u mt are divergence free and that u BL is located 
in a boundary layer of size 5e (where E stands for 'Ekman') near the 
surface. Denote by A BL ,A BL the size ofu BL , u BL in L°° , and by A int 
the size of u mt in L 2 (uj). 
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The assumptions above entail that A BL = A int ; on the other hand, 
since u BL is divergence free, we have 

aBL _ 1 aBL _ ^ Aint 
h — ~T~ 3 _ ~7~ A ■ 
OE Oe 

Consequently, since 

u BL ~A BL exp 

we infer that 



5 E 



i Bin IT~ aBL 1 Aint 

\ u h IIl 2 (u,) = V°EA h = ~jj^ A ■ 



Thus the energy in the boundary layer is always larger than the energy 
in the interior of the fluid with this type of model. The assumption 
that ||ug^* =1 || ~ ll w ftja=i[l s ^ ems from observations of the isothermal 
surfaces in the ocean, as we will explain in the next paragraph. From 
a physical point of view, having \\u BL \\ much larger than ||u m '|| is in 
fact quite reasonable: indeed, it is observed that subsurface currents 
generally travel at a much slower speed when compared to surface flows. 

Let us also emphasize that in the case of the /-plane model (i.e. when the 
rotation vector b is constant), the result of Theorem 11.11 is false in general. 
Indeed, the interior part of the solution must satisfy the geostrophic system, 
namely 

+ V hP = 0, 
div h u h + d z u 3 = 0, 
d zP = 0, 

and thus u is a two-dimensional divergence free vector field. In other words, 
■U3 = and thus the Ekman pumping velocities must be zero at first order. 
Consequently, the interior part of the solution cannot be wind-driven at first 
order. 

We now address the question of the stability of the stationary solution 
constructed above: 

Theorem 1.3 (Waves associated with equation (jl.7p ). Assume that ujh = 
TxR, and that b{xh) = 0y for all Xh = (x, y) £ uj^- 

For any e > 0, let v e be a solution to the propagation equation 

d t u + -b(x h ) Au+ ( 7q P J - v h /\ h u - ed zz u = 0, 

(x h ,z) Gu h x (0,1), 
with Vh = 0(e 2 ), supplemented with homogeneous boundary conditions 

d z u h \ z=0 = d z u h \ z=1 = 0, u 3 \ z=Q = u 3 \ z=1 = 0. 
Then v e can be decomposed as the sum of 
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• a stationary part v e (t, y) = J T Jq v e (t, x, y, z) dx dz, which satisfies 

8 t v e - v h d 2 y v e = 0, 

• Rossby waves v R = J v e dx% — v e corresponding to the 2D vorticity 
propagation 

d t c R + ^d x Aj; 1 e R -v h A h ( R = Q, 

where Cr = rot^ v R , 

• and gravity waves v e G = v e — J v e dx^. 

Rossby and Gravity waves have a dispersive behaviour as e vanishes: 

• Rossby waves disperse on a small time scale 

Vi>0, \/K CC lo, \\v R (t)\\ L 2 {K) -> as e ->■ 0, 

since we have assumed that y G R; 

• Gravity waves generate fast oscillations with respect to y, which slows 
down the propagation 

Vif CC uj, \\v e G (t)\\ L 2 {K) ^0 as (e,t) -» (0,oo). 
Remark 1.4. 

• Notice that the energy associated with gravity (or Poincare) waves 
propagates on a time scale much larger than the one of Rossby waves. 
This is due to the thin layer effect, which causes the apparition of 
small scales in the variable y. 

• The field v e is said to be "stationary" because the horizontal viscosity 
Vh is small: hence 

v e (t) f« vf t=0 in L 2 

on time scales of order one. 

Corollary 1.5. Assume that uJh = T x R, and that b(xh) = Py for all 
Xh = (x,y) G cuh- Assume that = 0(e 2 ). 

For any e > 0, let u e be a solution of (jl.ip supplemented with (jl.2p - (jl.3p . 
and assume that 

sup ||<| t=0 - <°*|U2 H + e||u|| t=0 - uf at \\ L 2 {ul) < +oo. 

£>0 ' ' 

Then for any finite time t > 0, 

u\t)-u stat ~v\t) + vh{t) m L 2 loc {u) 

where Vq is the (slow propagating and fast oscillating) gravity part of the 
velocity field v e defined in Theorem \1.3\ and v e is the stationary part of v e . 

(Note in particular that the vertical component of the velocity u t is not 
expected to be bounded - as is usually claimed for shallow water approxima- 
tion.) 

The above Corollary is an immediate consequence of Theorems 11.11 and 
11.31 together with the energy inequality. 
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1.3. Towards a mathematical derivation of the thermocline. 

In this paragraph, we try to justify the shape of the surfaces of equal 
temperature in the ocean, in view of the results of Theorem 11.11 

The isothermal surface which is located just below the Ekman bound- 
ary layer is of special interest to oceanographers, due to its importance on 
the global oceanic circulation (see [23l [24J, [18] ) . Figure CD below shows the 
longitudinal variations of the temperature in the Pacific ocean in a layer of 
1000 m depth below the surface. In particular, there are zones in which 
the temperature surfaces are tilted up (that is, there is a flux of cold water 
towards the surface); this phenomenon cannot always be accounted for by 
the heating differences at the surface, as shows the upward flux of cold water 
in the equatorial zone. The physical justification of these particular shapes 
is the following: inside the ocean, the temperature T solves an equation of 
the kind 

u* ■ VT - kAT = 0, 

where u* is the velocity of oceanic currents (in the dimensional variables) 
and k is the heat conductivity coefficient. If the temperature diffusion can 
be neglected, this equation takes the form 

u * . VT = 0, 

which means that u is a tangent vector to the isothermal surfaces. Con- 
sequently, the temperature surfaces are tilted up (or down) if and only if 
u 3| SU rface ^ °> or more precisely, if \u 3 \ z=l \/\u h \ z=1 \ = O(l) in rescaled vari- 
ables. This justifies the assumption 

ttn.int II Hn.in-t II 
|P3|*=lll ~ ll u h|«=lll 

in the previous paragraph (see Remark 11.21 (iv)). 

In that regard, the special solution constructed in Theorem II .11 is of par- 
ticular interest. Indeed, in rescaled variables, we have (see section [3]) 

%| 2 =i = -^rot ft (7-^<7i. 

Hence U3| sur face 7^ 0, and our model predicts that the temperature surfaces 
are indeed modified by the Ekman pumping velocity. 

We now give a rigorous result about the asymptotic shape of the tem- 
perature in our model. We denote with a star the original variables. We 
write 

T(xl,z*)=T + T 1 e(^, Z - 

with the same notations as in paragraph 11.11 The temperature To is a 
reference temperature (for instance, To = 10°C), whereas T\ is the order 
of magnitude of the variations of the temperature. Performing the same 
change of variables as in paragraph II .1^ we obtain 

u ■ V0 - \r?A h 9 - Xd zz 9 = 0, 
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Figure 1 . Longitudinal section of the surfaces of equal tem- 
perature in the Pacific ocean (from the WOCE Pacific Ocean 
Atlas). 



where the diffusion coefficient A is given by 



We recall that r/ is the aspect ratio of the domain; as in Theorems 11.11 and 
11.31 we take rj = e. With the notation of Theorem 11.11 our result is the 
following: 
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Proposition 1.6. Let A > 0. Assume that the wind stress a S H s (uh) is 
such that 

\a(x,y)\ <Cy k V(x, y) G u h , 
\Va(x,y)\ < C\y\ k ~ l V(x,y)eu h , 
for some k,s>2 chosen sufficiently large, and assume that 

(1-14) [|V fc «^|U<-( w) < - 

Let 9 be the solution of the equation 

(1.15) u stat ■ V9 - n 2 XA h 9 - \d zz 6 = 0, 
supplemented with the boundary conditions 

(1.16) 0\z=i = Oi, d z 9 {z=0 = 0, 

for some function 9\ E H 2 {u>h). 
Define the function 9 app by 

9 a PP(x h , z) = 9(x h , z) + e9 BL (x h , ^ 

where 9, 9 BL are solutions of 

-\d zz 9 + u mt • V0 = inu, 

(1.17) ! 

9\ z =i = #i, d z 9\ z=0 = 0, 

and 

-Xdcc9 BL {x h , C) + eu BL (x h , 1 - eC) • V/A = 0, 



■\BL 



(x h ,Q — ► o. 

f — »oo 



Then as e — ► 0, 

ll«-e <w ||^ (w) + ||5 z ((9-^)|| £aH -.0. 



Remark 1.7. (i) T/ie assumption (|1.14|) on i/ie size o/ V^u^ 1 ' is purely 
technical, and does not have any physical interpretation. It rises from 
the fact that equation (|1.17|) on 9 is degenerate in the horizontal vari- 
able; we refer to section^ for more details. We emphasize in particu- 
lar that if (|1.14[) is not satisfied, equation (|1.17|) is still well-posed in 
L 2 (ujh,H l {Q, 1)); however, in this case, we are no longer able to prove 
the convergence. 

(ii) The boundary conditions (I1.16P mean that the atmosphere acts like a 
thermostat for the ocean, and that there is no heat flux at the bottom 
of the ocean. Both assumptions seem reasonable from a physical point 
of view, although other boundary conditions might also make sense: 
for instance, it could also be assumed that the heat flux at the surface 
is a given function of the latitude. 
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(iii) Let us mention a last direction towards which the physical accuracy of 
our model could be improved. When considering the spatial variations 
of the temperature, it would be more reasonable to consider a model 
which couples the velocity of ocean currents and the temperature, in 
the spirit of [3j. However, the relevant scalings within such models 
are not completely clear. Furthermore, the analysis in Chapter 6 of 
[23j shows that for such problems, the curvature of the Earth should 
be taken into account. Hence we leave this issue aside in the present 
paper. 

The construction of the article is as follows: in the next two sections, 
we construct the stationary solution of equation (jl.ip . starting with the 
boundary layer part, and then building the interior part by solving the 
geostrophic equations with a Dirichlet boundary condition on the vertical 
component. Then, we prove Theorem 11.31 in sections U] and [5l by treating 
separately the two-dimensional and three-dimensional parts of the initial 
data. Eventually, section [6] is dedicated to the proof of Proposition 11.61 

2. The boundary layer part of the stationary solution 

In this section, we construct functions u BL ,p BL which are approximate 
stationary solutions of equation (II. lh (in the sense of Theorem II. ip , and 
which satisfy the horizontal part of the boundary condition (|1.3p . These 
functions are located in a boundary layer in the vicinity of the surface z = 1. 
Our methodology is the following: we first assume that vh = 0, and we use 
the classical construction of Ekman layers in this case. We then derive 
several estimates on the functions thus obtained. Eventually, we estimate 
the error terms in equation (II. ip which are due to the fact that vh is non 
zero. 

2.1. Construction in the case = 0. 

When the horizontal viscosity vanishes, the construction of the boundary 
layer is exactly the same as in the /-plane model, i.e. when the function 
b does not depend on Xh- Indeed, in this case the variable Xh is merely a 
parameter of the equation, and building the boundary layer term amounts 
to solving an equation on the rate of exponential decay. For more results re- 
garding classical boundary layers, we refer to [H 120 } I2T1 125] . Nonetheless, let 
us stress that even though the construction itself is the same, the estimates 
become much more involved than in the case of the /-plane model. Indeed, 
the vanishing points of b create singularities, and prevent the boundary layer 
terms to be in L 2 in general. Hence, assumptions on the stress a have to be 
introduced in order to handle these singularities. 

The construction of the boundary layer term is as follows: we wish to 
construct an approximate solution (u BL ,p BL ) of (II. lh . such that (jl.3p is 
satisfied. Furthermore, we assume that this approximate solution is small 
outside a boundary layer located in the vicinity of the surface z = 1. Hence, 
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,BL „BL 



we look for u ,p in the form 



»// 7 - (/..'■/,.:)■ ■ If 1 {■<hX-^- 

1 - z 



p BL (t,x h ,z)=P BL ( x h , 



We assume that U BL , P BL together with all their derivatives vanish as £ — > 

00, where £ stands for the rescaled variable (1 — z)/e. Inserting the above 
Ansatz into equation (jl.ip yields 

( b(x h ){U BL )^-d}U BL + eV h P BL = 0, 

(2.1) J -^<7 3 Bi - ^d c P BL = 0, 
[ div ft l/f L - ±d c /7 3 = 0. 

The last two equations entail that 

pBL = _ e 2g cU BL = _ £ 3 ^ jjBL 

We henceforth neglect the pressure term in the equation on U BL . Then, we 
set, as usual (see for instance [20]). 

U£:=U h ±iU£. 

Above and in the rest of the article, for all u = (1*1,1*2) £ R 2 5 u 1 - := 
(-u 2 ,ui). 

An easy calculation leads to 

Consequently, is an exponentially decaying function of the form 

u h(x h ,0 = -r+7 — ±icr ± )(x h )exp(-\ ± (x h )(), 
eA {Xh) 

where the decay rate X^ 1 is defined by 

(A±) 2 = Tib and 3?(A ± ) > 0, 

1. e. 

(2.2) xH*>) = x t M = 1 * i % m m?»- 

Notice in particular that the decay rates X^ 1 vanish at y = and depend 
only on y. 

Going back to the definition of £7. , we infer that 

(2.3) !/«(*», = !JH = ij; (fL±igM eX p(-A*(x k) c). 
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Hence, in order that U BL is divergence free, we set 

/DC 
div h U* L (xh,C')dt' 

~\ E( div ^ =F "ot^)(x ft )(A ± (x ?t ))" 2 e - A± ^) c 



2 ± 



(2-4) + \ ^> ± ia X ){x h ) • (2 + CA ± (x fc )) 



)e - A± (^)^. 



± 

We have used the convention 

rot^ u h = - div h 

for two dimensional- vector fields. 

The remaining flux term is then given by 

(2.5)u^ =1 = U.^ =Q (x h ) = -lYs^haTirothaKx^ixh))- 2 

We now wish to point out a particular difficulty steming from the above 
construction. If the Coriolis factor b has vanishing points, which occurs 
in particular in the case of the /3-plane approximation (b(xh) = fly), then 
the functions U B L , U BL may not be square integrable if the function a is 
arbitrary. Hence, the function a should vanish at a sufficiently high order 
near y = so that the singularity disappears. We will check that (II. 8p entails 
that the functions U BL ,U BL defined by (|2.3p . (|2.4p are square integrable. 
For further purposes, we also require that the function V 'hU BL belongs to 
L 2 (tOh x [0, oo)^). Unfortunately, assumption (|1.8p is not sufficient to ensure 
such a result. Thus we introduce an approximate boundary layer term, in 
which the low values of b have been truncated. 

2.2. Estimates on the boundary layer terms. 

We begin with a short justification of the need for a truncation. Using 
the definition of \^ together with assumption (ILlOh . we infer that if y is 
close to zero, then 

IIV7 TT BL (n- \\\ < rf HZMgftll | \ Va ( x h)\ | W(Xh)\ \ 

\\vx h u 3 (^)|| L 2 ([0)Oo)f) < c^ y5/4 + ^ 9/4 

Hence X7 Xh U BL does not belong to L 2 (a^ x [0, oo)^) in general. We thus 
define, for any 5 > 0, the function 

(2-7) b s (y) = b(y)if>0^\, y^O 
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where ip S C°°((0, oo)) is such that 

ip{y) > 2 for yG (0,oo), 

iP(y) = l if 2/ > 2, 

V(y)=y- Q if ye (0,1), 

for some exponent a £ (0, 1) to be chosen later on. Notice that with this 
choice of ifi, the function b$ behaves like 5 a y 1 ~ a for y > near zero. Conse- 
quently, bs vanishes with a weaker rate than b, and thus cr/bs vanishes more 
strongly than a/b. 

We now define approximated decay rates Xf by replacing b by b$ in the 
expression (12. 2ft ; eventually, we define approximated boundary layer terms 
by the formulas (|2.3j) - (12.4p . in which the decay rates X^ 1 have been replaced 



by A± 

We then have the following result: 

Lemma 2.1. Assume that hypotheses (|1.8p . (jl.lOj) are satisfied. Then there 
exists a constant C , depending only on a and b, such that for alia > 0,5 > 0, 

\\ttBL\\ ^ C 

\\ u S,h IIZ2(w fc x[0,aa) c ) S — 

1 1 C7 ^3 Z ' 1 U 2 X [0, 00) c ) < C- 

Additionally, if a > 3/5, there exists a constant C a , depending only on 
a, a and b, such that for all 5 > 0, 

\\X7 tt bl \\ <r C° 

\\VhU5,h ll£ 2 (u> h x[0,oo) c ) S — 

WhU 5 ,3 ||L2( Wfc x[0,oo) c ) < -p/l- 

Moreover, for all 5 > 0, 

£11/4 

- 65)^,/! llL 2 (a; h x[0,oo) f ) < C— — , 
llVfePj ||L2(a;, i x[0,oo) f ) < Ca^yJ- 

Remark 2.2. TTie above estimates are given for the rescaled boundary layer 
profiles U BL , P BL , which are defined on x [0, 00)^. Remember that the 
boundary layer part of the stationary solution is defined on lu^ x [0, 1] by 

u BL (x h ,z) = U^ L (x h ,- 



Hence 

\\ u h llL 2 Kx(0,l)) - 6 \\ U h llL 2 (^x[0,oo))- 

The same estimates hold for p BL ,u B . 
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Proof. • L 2 estimates: According to (jl.lOp . and to the definition of bs, we 
have 

|A ± (x ft )| /Oify/0, |A ± (x ft )| ~ T^jasy^O, 
and thus there exists a constant C such that 

|»(A ± (or h ))r 1 , |A ± (x h )|- 1 <C| I /|- 1 /2 Vx, Vy G [-1,1]. 
Similarly, for all 5 > 0, we have, for |y| < 5 

If 1 2/ 1 > <5 5 then Ac (x^) satisfies the same estimates as A ± (x/j). A careful 
computation leads to 

(2.8) y o TO^oi ^ = ^k(*h)i 2 x;ui — : — 



± 



A±(y)|^(A±(y))- 



Hence we obtain 



{l i 5+3" t— 2^ . i i 

|y| 4 S 4 if 1 2/| < 5, 

M 5/4 if<5<|y|<l, 

Eventually, we infer that 



la I else. 



\TT BL \\ S 
\U6,h IU 2 Kx[0,oo) f ) S ( 



where the constant Co depends only on b and a. Notice that the truncation 

J B 



does not play any role at this stage: the same arguments show that U? L G 



L 2 (uj h x [0,oo) c ). 
Similarly, we have 



( 



2.9) r\ Us f M \ 2 d ( <cj:^L + 

Jo V A a (y)| 5 



Va(x h )\ 2 , \a(x h )\ 2 \V\f(y)\ 



2 



\K(y)\ 7 



Using the definition of the decay rates together with the definition of the 
function ip, we obtain 



^5 a / 2 if \y\ < 5, 




else. 



Thus 

T BL 



U 8,3 llL 2 (o; h x[0,oo) c ) < Co- 



i?^ estimates: 



We begin with the bound on V^t/pif; the calculations are very similar to 
the ones which led to the 1? bound on Us 3, and are therefore left to the 
reader. In fact, the situation is even a little less singular than in the case of 
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Ug3 (we "gain" one integration with respect to the variable £, and thus one 
factor (A^)" 1 ). The bounds on A^ 1 and a entail that 

\\X7 TjBL\\ ^ Co 

\\Vh U 5,h llL 2 (o; h x[0,oo) c ) S — • 

We now tackle the bound on VhJJf^: first, differentiating equation (|2.4p 
with respect to Xh, we obtain 

IIV^WIMo,^ < C (^ + |V, W |Jg«) 

In the expression above, we have denoted by |A<j| the common size of |A^| 
and |A7|. Notice that due to the sign change in b at y = 0, there is in general 
a Dirac mass at y = in the term d yy \s; more precisely, the part of V^J/^g 
which is not absolutely continuous with respect to the Lebesgue measure is 
of the type 

*v=oMjb^p = ^=oklN- 9/4 - 

At this stage, the need for a truncation is clear: if b$ is replaced by b, then 
|cr||6| _9//4 ~ |y| _1//4 near y = 0, and thus the singular part of VhU BL is not 
well-defined in the sense of distributions. Conversely, if a > 1/9, then 



9a-l „ 9a 



\cr\\bs\ ' ~ |y| 4 6 4 as y — > 

and thus the singular part of ^hUf^ is zero. 
Gathering all the terms, we deduce that 



4^4 if |y| < d, 

l|V h C^(x fc )|| L a<C^ | y |-5/4 if<5<H<l, 

|cr(a;/,)| + \Va(x h )\ + \D 2 a{x h )\ else. 



Thus ^hUfs £ L 2 (ujh x [0, oo)) if and only if a > 3/5, and in this case there 
exists a constant C a , depending on a, b and a, such that for all 5 > 

C 

• Error estimates: First, by definition of we have 
\\(b-b s )U s f\\l 2 = / r\b(y)-b 5 (y)\ 2 \U s B /t(x,y,()\ 2 d(dydx. 

Ju h n{\y\<28} JO 
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Notice that for all y 6 R \ {0}, 
\b(y)-b s (y)\ = \b(y)\ 
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l\y\<s\b(y)\ 



y\ a 

< cii^slyl 1 -"^ 

1— a 



1 ] + l 5 <| y |<25|6(y)| 

" |yD + 1 5<|y|<25|y|) 



< Cl|y|<2<5|y 

Using (|2.8p . we infer 



< 



< s 

e z 



\y\ 2(1 - a) S 2a \y\'\y\ 

x£T J\y\<26 

y \(9- a )/2 6 a/2 dy 



-3/2 



3a/2 



< 



y\<2S 
C5 11 / 2 



There remains to evaluate PP L ■ By definition, 



BL 



-e 3 div h U, 



BL 



Using the same kinds of calculations as the ones which led to the bound on 
VhUf%, we deduce that 

■2 



\^h P 5 llL 2 (u; h x[0,oo) f ) < C 



s 1 / 4 ' 



□ 



2.3. Error estimates in the case 1/^^=0 and conditions on the pa- 
rameter 5. 

If Vh 7^ 0) we keep the construction of the previous paragraph, and we 
merely treat the viscous terms as error terms. The function uffc is an 
approximate solution of the horizontal part of equation (jl.ip . with the error 
term 

bsXuft) 1 - - v h A h uf£ + V h pf L . 

According to the estimates of the previous paragraph (see Lemma I2.1|) , we 
have 



< C 



L*(u>) 



e 3/2 



and 



\WhA h U Sth ||L2([0,l],H-i( Wh )) < C — 

,5/2 

„BL\\ _ / n e 

<5V4" 



\\v hP ^\\ LHul) < c 
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Recall that because of the boundary layer scaling, there is a factor e 1 / 2 
between the L 2 norms of uf and Uf L . 

Hence, in order that the conditions of Theorem 11.11 are satisfied, the 
numbers e,Uh,5 should verify 

(2.10) e 10 «5«e 6 / n , ^«e. 

On the other hand, uf L is an approximate solution of the vertical com- 
ponent of equation (jl.ip , with an error term equal to 

-u h A h u BL . 

The estimates of the previous paragraph entail that 

\\-v h A h ui L \\ L2mhH „ 1(ujh)) < C-^-. 

In order that the condition of Theorem 1 1.1 1 is satisfied, the parameter 5 must 
be chosen so that 

(2.11) 6 » v 2 h /3 e 2 . 

Notice that if = O(l) and e = o(l), we always have 

^ / V«6 6 / 11 . 

Hence it is always possible to choose a parameter 5 which matches the above 
conditions. 

Further conditions on the parameter 5 will be given in the next section. 
When the conditions (|2.10p . (|2.1ip are satisfied, the couple {uf L ,pf L ) is an 
approximate solution of equation (jl.ip in the sense of Theorem 11.11 Fur- 
thermore, uf L satisfies the horizontal part of the boundary condition (|1.3p 
at z = 1; u BL , on the other hand, does not satisfy the non penetration con- 
dition at z = 1. Hence, we construct in the next section an interior term, 
which is also an approximate solution of (jl.ip . and which lifts the trace of 
u BL at z = 1. 

Notice that u BL also has a non-vanishing trace at z = 0; however, this 
trace is exponentially small on the set where b is bounded away from zero, 
and can thus be lifted thanks to an exponentially small corrector. This will 
be taken care of after the construction of the interior term u mt , in the last 
paragraph of the next section. 

3. The interior part of the stationary solution 

In this section, we construct a stationary solution u mt of equation (jl.ip . 
which is such that u mt + u BL satisfies the boundary conditions (jl.2l) . (jl.3p . 
Going back to equation (jl.ip . it can be readily checked that the function 
u mt should satisfy the system 

%)«) ± + V^ = 0, 

(3.1) d z p = 0, 

divu int = 0, 
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together with the boundary conditions 



.int n „.int „.BL 

(3.2) 



a ,.%nt _ n int _ l 



®z U K\z=Q — °> U 3\z=0 — °- 

We recall that since the function vM L depends on the small parameter S, the 



function u m also depends on 5 in general, and thus will be denoted by uf lt 
in the sequel. Hence we also investigate the asymptotic behaviour of u 1 ^ as 
5^0. 

It turns out that the solution of the sytem (|3.1jl - (j3.2j) is unique, up to a 
function of the type (v(y),0, 0). Hence we give in this paragraph a straight- 
forward way of building the solution, and then we derive L 2 estimates on 

in 
5 



the function uj lt . The main result of this section is the following: 



Lemma 1. Assume that assumptions (ll.8p -( fTTT0l) are fulfilled. Then there 
exists a solution uj 11 G L 2 (uj) of the system (13. lj) . Moreover, there exists a 
positive constant C , depending only on a and b, such that 



\K lt \\ L ^)<c v^>o. 

3.1. Construction of u l s nt . 

To begin with, we differentiate the first equation of (13. ip with respect to 
z, and we obtain 

b(y)dM n h) L = 0- 

Since u 1 ^ 1 is divergence-free, we infer that d zz ii!fi = 0. Hence the third 



.int 
l S,3 

,BL 



component u^S is uniquely determined; in order to lighten the notation, set 



w s (x h ) = -u S3lz=1 (x h ). 
We have 

Us$(x h ,z) = zw 5 (x h ). 

Then, taking the two-dimensional curl of the first equation in (|3.ip . we 
derive 

rotfcCftCug)- 1 ) = dhr fc (6«g) = 0. 
Since the Coriolis factor only depends on the latitude y, we are led to 

b'(y)utS = -6(y)div h ug = +b( V )d g v$ = b(y)w s (x h ). 

Consequently, the second component is also uniquely determined. In the 
case when b{y) = f3y, one has in particular 

*4$(sh) = yws{x h ). 

This equation is known as the Sverdrup relation (see \23\ 124] ). 
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There remains to compute the first component of u mt ; the divergence- free 
condition entails that 

d x V$ = -dyU^i - d zU fi = -d y (yW S ^j - W 5 



<>» (£) 



9 bb "\ b X 

Wj ws ~ V yW5 ' 



Notice that this equation has a solution in if and only if the right-hand 
side has zero average in x, for all y. This is satisfied in particular if 



(3.3) / w s (x,y) 

Jt 



dx = Vy. 



We assume that this assumption is satisfied for the time being, and we will 
prove that it is in fact equivalent to (|1,9|) . Integrating the equality giving 
fljUrj with respect to x, we deduce that u\ n ± is defined up to a function of 
y only, provided (|3.3p is satisfied. 

Now, let us compute ws in terms of a and b. Using equation (|2.5p . we 
infer that 



ii'M .)•/, i - - ^ (div h a =F i rot h a) , - ^(cr ± ia~ 



2^— — '(Aft. ' (A±)3- 

By definition of A 1 * 1 (see ()2.2p ). we have 

± _ / 1 =F » sign(fr) sign(fe)^ 

Hence 

(3.4) Ws (x h ) = I rot, a + 1^ • V6 5 = ^ - d y 

(Recall that 65 only depends on the latitude y.) 

We now prove the equivalence of (jl.9p and (|3.3p . It is clear that (jl.9p =^ 
(|3.3p . Conversely, if (|3.3p is satisfied, then (|3.4p leads to the existence of a 
constant as £ R such that 

Since o\ vanishes quadratically near y = 0, we deduce that the left-hand side 
of the above equality vanishes at least linearly near y = 0. Consequently, 
as = for all 5, and thus (II. 9p is satisfied. 
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3.2. Bounds on u int . 

We begin with a bound on the function w$ given by ()3.4j) . We recall that 

~ 0v near y = °> 

and 

ax(x,y) = -S^a 1 (x,0)y 2 + O(y 3 ) as y -► 0, 
d x a 2 {x,y) = 0(\y\ 2 ) as y -> 0. 

Thus 

/(Jl\ b$dy(Tl ~ Oldybs . n I I o 

j v 65 J = — — P — = ) y ~* ' - 



The exponent a was introduced in the previous section, see (|2.7p . 

Consequently, there exists a constant C (independent of (5) such that 

IKIIlv/J ^ a 

This entails immediately that u s n ^ and are bounded in L 2 (uj), uniformly 
in 5. 

As for uf 1 *, we have, by definition 

a int a fbd x a 2 \ d x a 2 , „ (b o\\ , „ o"i 



4 s — u^-a 



Integrating with respect to the variable 2, we deduce that 



<--*UFj-«7 + *#-*U(wi-\ 

where Si(x,y) = a±(x',y) dx' . Using the definition of bs, we obtain 



b% sv^y 

It can be checked that the function in the right-hand side is C°° on (0, 00) 
and bounded, together with all its derivatives. Moreover, its support is 
included in [0,25]. As a consequence, the term 

a, ill (5-1 



b s \b'b s 

is o(l) in H 1 ^) as 5 — > 0. The other terms can be evaluated in a similar 
fashion. Using the assumptions on a and b together with the definition of 
"0, we deduce that there exists a constant C[a] such that 

(3-5) IK n 'lU» < C[a]. 
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• We now derive estimates in L 2 ([0, 1], H (a^)), which are needed to 
bound the error term UhAhuJ 1 *. First, using the definition of bs together 
with assumptions (11. 8ft . (jl.lOj) . it can be proved that 

d y w 5 = 0{y a -H- a ) as y 0, |y| < 6. 

Hence d y ws G L 2 (lo) (recall that a>3/5>l/2) and 

\\d y ws\\ L 2 (Lu) = 0(<T 1/2 ). 

The term d x w$, on the other hand, is bounded in L 2 (lo), uniformly in 5. 
Consequently, there exists a constant C, depending only on a, b and a, such 
that 

C 

\\Vhv%£\\ L 2 {uj ) < ^r/2- 

Similarly, we prove that d y us : 2 = 0(\y\ a 5~~ a ) for y in a neighbourhood of 
zero, and thus there exists a constant C such that 

IIV^II^H < C. 

We now tackle the term ug i) using either the expression of d x us t i in terms 
of ws or the final definition in terms of o<i and Si, it can be checked that 

d y u 5> i = 0{y a - l 5- a ) as y - 0, |y| < <J. 

The largest terms are those coming from S\ (or from bd y ws/b'); for instance, 
the above calculations show that 

ydyw 6 = 0(\y\ a 5- a ); 

since one power of y is lost with each differentiation with respect to y, we 
obtain the desired bound on ug^. Eventually, we are led to 

C 

• Notice that 

as 5 — > 0, where u mt is the function defined by the same expressions as uj 11 , 
but replacing every occurrence of w$ by 

rot/i a a 1 - • V6 

By definition of bs, w and ^5 coincide on the set {|y| > 25}. Moreover, w is 
bounded in 1? and w, yd y w have finite limits as y — > 0, while 



\ w s\ + \y\ \d y w s \ = 0(1). 

\y\<S 

Consequently, ws (resp. bd y ws) converges towards w (resp bd y w) in L 2 (u>) 
as (5 — ► 0. The convergence of u^ 1 ' follows. However, in general, u int does 
not belong to H 1 {oj), except if the surface stress a vanishes at sufficiently 
high order. 
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3.3. Proof of Theorem QUI 

Let us first evaluate the error terms in equation (II. 7ft . To begin with, 
notice that d zz vH'^' t = 0, so that there is no error term associated with the 
vertical Laplacian. Consequently, the only error terms in equation (II. lj) are 
those coming from the term i^A/jit^ ■ . 

According to the 1 estimates of the previous paragraph, we have 

II 7/ A 7,* n *ll <r n Vh 

Wh^hU^ h \\ L2miH _ 1{uih)) < G-^=, 

II 7/ A ,,* n * II <r n Vh 

\\^hU 5 , 3 |y ([0)1 ] ifr -i W) -°7|- 

In order that the conditions of Theorem 1 1.1 1 are satisfied, we have to choose 
the parameter (5 so that i > i/jj. We recall that (5, should also satisfy 
(|2.10p . (|2.11|) . Thus the new conditions on 5, are 

(3.6) max(^,e 10 ,^ /3 e 2 ) <5<e 6/n , ^ < e. 

• The proof of Theorem II .11 is now almost complete. There only remains 
to take care of the boundary conditions: indeed, as we have explained at 
the end of the previous section, the trace of d z uffc and uf^ is non zero at 
z = 0. Hence, we define a corrector v$ nt , which is small in H , and which 
lifts the remaining boundary conditions. The result is the following: 

Lemma 3.1. Assume that Vh = o(e) and that there exists k £ (1,2) such 
that 5>e K . 

Then there exists a divergence free function vf 1 * ' , such that vf 1 = o(l) in 
L 2 (lo), which satisfies the conditions 

d*v$ {z=1 = 0, v^ z=1 = 0, 

°z V 5,h\z=0 ~ ° zU h\z=0i v S,3\z=0 ~ u 3\z=0- 

Furthermore, we can choose the parameter a of the truncation function tj) 
so that 

he 3 A vf\ -ed zz v{ nt = o(l) in L 2 {uj), 

-ed gt vt a = o(e- 1 ) inL 2 (u;), 
- h ^ h vl nt = o(l) m L 2 ([0, 1]), H-^Uh)), 
- h ^ h vf l = o{e~ l ) m ^([OdD.ff^K)). 

Before proving the lemma, let us complete the proof of Theorem 11.11 we 
choose a parameter 5 which matches the conditions of Lemma 13.11 together 
with (|3.6p . Notice that the choice 5 = e works. We set 

..stat _ BL i int , int. 

U —Us t U S "+~ v 5 > 

by construction, u stat satisfies the boundary conditions (|1.2|) . (|1.3|) . and it 
is an approximate solution of equation (jl.ip . The bounds on u™ 1 and uf 
were proved in the previous paragraphs. Hence Theorem 11.11 is proved. 
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Proof of Lemma \3.1[ Throughout the proof, we drop all indices 5 in order 
not to burden the notation. 

The construction of the corrector v mt follows the one given in Lemma 1 
in Appendix B of [9]: setting 

<f>h := -d z u^ L z=0 , 3 := -uf z=Q , 

we define 

v t Z* = ^- r L <i>h + Vhx, 

where the potential x £ B ( u h) 1S defined by 
AftX = J div h vf l - - div h <j) h 

= - [vi nt ] z z zl ~ \ div h <j> h = 3 - ~ div ft (j) h . 

o b 

We will check later on that the function 03 has zero mean value on Uh, so 
that x is well-defined. The third component of v int is then determined by 

vf\x h , z) = - f d z vf\x h , z')dz' = f div h v{ nt (x h , z') dz 1 . 

J z J z 

By construction, v mt is divergence free and satisfies the correct boundary 
conditions. There remains to evaluate v mt in L 2 (u>) and L 2 ([0, 1], H 1 (uJh))- 
The boundary conditions 0^, 03 are given by 

A±' 



2e^ V 

(.).-.. -- ^(div ft n + i rot/, T)-7-r- r exp 



2e^- - m J6XP V e 

1 / A± 

2^ (A*) 2 V e 



± 

Recall that in the expressions above, the functions A 1 * 1 are in fact A^. Notice 
that 

03 = div h ip, 

where 

1^, I 1 / X ± \ 

l P = -2^ (T±ia )(A±)2 exp V" —) ; 

this proves that 03 has zero mean value on uJh, and will be used several times 
in the proof. 

We now derive three type of estimates: first, estimates of div/j 0^ and 03 in 
L 2 {uih) will yield i/ 2 (u;/ 1 )-bounds on x, and thus bounds in L 2 ([0, 1], H l {u}h)) 
for the function ut ri *, and in L 2 {ui) for the function ul^ 1 *. Then, estimates 
of 4>h and f will provide L 2 (u;)-bounds on uj^*. Eventually, L 2 estimates of 
\7ft03, -D^0/i will allow us to derive bounds on v| nt in L 2 ([0, 1], iT 1 ^)). 
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• Estimates of divh^h and ^3 in L 2 (uJh)'- 

The main difficulty lies in the fact that Xf does not have the same be- 
haviour for \y\ < 5 and \y\ > 5. We merely explain how the term div^ (p^ is 
evaluated; the treatment of the term ^3 is left to the reader. 

If |y| > 1 5 we have 



(A 



±\2 



Tib(y), with \b(y)\ > C (see (fTTTOl) h 



Thus 



and 



exp 



A± 
e 



< exp 



C 
e 



[ [ \ div h c/) h (x, y)\ 2 dxdy <^\\a\\ 2 m( ) exp (-—) . 

J\y\>lJxeT e y ' \ e J 

On the set where 8 < \y\ < 1, the assumptions on the truncation function 
■0 entail that there exists a constant c such that 

c" 1 !?/! 1 / 2 < 3^(A±(2/)), |A±(2/)| < c^l 1 / 2 , 
\d y \f(y)\ < c\y\- 1 / 2 . 



As a consequence, 

I 



C 

* T 2 



&<\y\<l JxGT 
1 



div h 4> h (x,y)\ 2 dx dy 



\y\ exp 



y rr ex P 



-2c^f) (hi 



There remains to treat the set where \y\ < 5; because of the truncation 
function tp, this part is the most complicated. The definition of the function 
ip and the fact that b(y) ~ (3y for y close to zero entail that 

eT^l^f < \Xf{y)\,U(Xf(y)) < c\y\^Sf , 



\d y Xf(y)\<c\y\ 



l + ct a 
2 2. 



Thus, for instance 







/ / 


div/jfiexp 


)\y\<8 JxeT 





c 



< c 



W\ exp 

2 \ 3 



\y\ 2 02 
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The other terms in div^ fa are evaluated in the same way. Gathering all the 
terms, we infer that 



divfcfo||i2 (a)fc) <C I e- 2 exp(-C/e) + e 2 + e 1 



e 



1— a 



The corresponding term in ^/uhAhV mt should be o(l) in L 2 ([0, 1], H 1 (ujf l )); 
hence the parameters e,Vh,5 must satisfy 

^(e- 2 exp(-C/e)+6 2 ) =o(l), 

V^e- 1 (-^-\ = o(l). 

It is obvious that for Uh, e <C 1, the first condition is always satisfied. The 
second condition reads 

Since j//,Ce (see (|3.6p ). we always have 

1 ~" 4+2q 5+a 

and 2 < (5 + a) /3a since a < 1. Consequently, provided that ^ <C e, we 
have 

1=2. 4+2c 9 

Hence, if S > e 2 , y^div/j^h = o(l) in L 2 (uJh)- It can be checked that </>3 
satisfies the same property. The same estimates also prove that 

ed zz vl nt = oier 1 ) inL 2 H, 

vi nt = o(l) inL 2 (u;). 

Similarly, we show that fa, cp = o(e) in L 2 (ujh) as long as 5 > e 2 , and thus 
v% h l = °( e ) m Notice that this is not entirely sufficient to prove the 
assertion of the Lemma if the Coriolis factor b is unbounded. However, 
using the fact that fa and <p decay like exp(— |6| 1//2 /e) for \y\ > 1, it can be 
easily proved that 

A h (b X ) = bA hX + 2b'd 2X + b" X = o(e) in H~ l (ui h ). 

Hence b\ = o(e) in H (u^), and bV\ = ^{^X)~^'x = °( e ) m Eventually, 
we infer that bv™ 1 = o(e) in L 2 (uj). 

• Estimates of .D 2 ^ and in L 2 (ujh): 

Calculations similar to the ones led above show that 

P^IU-w < c Mh , m {^lA + 1 + (^) *) • 

And if 8 > e 2 , -C e, then the right hand side is o(e^ 1 u h 1 ^ 2 ). 
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The term Vh4>3 is the most singular of all, and eventually prevents us 
from taking 5 > e 2 ; indeed, it can be proved that 

2(2q-l) \ 

e 2 5— J 

In order that the right hand side is o((e- v /P^) _1 ), we must have 

5a-l 

5 > e 2c . 

Since 5 > e K for some k 6 (1, 2), we choose a € (0, 1) such that 5c ^~ 1 > k. 
We then infer that Vh<t>z = o((e v ^) _1 ) in L 2 (w/ l ), and thus y^A/j^"* = 
oCe-^iii^aO,!],^- 1 ^)). 

□ 



\\Vhfa\\L2(u h ) < C| 



(7 



4. Two-dimensional propagation 

We recall that throughout this section and the following, we assume that 
b{xh) = (3y, and that u/, = TxR. The object of this section is to prove the 
"two-dimensional part" of Theorem 11.31 In particular, we prove that a two- 
dimensional perturbation of the solution u stat creates waves, propagating at 
a speed of order e _1 , with frequencies given by 

P k 

p \k\ 2 + \t v r 

where (k,£ y ) is the wavelength. 

A consequence of our result is that if u stat is initially perturbed by a 
two-dimensional function u° such that u° = 0(1) in L 2 and such that the 
x-average of u° is zero (i.e. vP has no Fourier mode corresponding to k = 0), 
then the solution of (jl.ip with initial data u stat + u° becomes close to u stat 
for finite times, with an error term which is o(l) in L 2 ([Tq,T] x uj) for all 
T > T > 0. 

Definition 4.1. Denote by P2D : L 2 (uJh) 2 — ► L 2 (uJh) 2 the projection on two- 
dimensional divergence free vector fields. The Rossby propagation operator, 
denoted by Lr, is defined by 

i/?^ = p 2D (6y ± ). 

Lemma 4.2. Let v9 £ L 2 (u>h) be a two-dimensional divergence free vector 
field, and let u £ C(R+, L 2 (uj)) be the solution of equation (jl.ip with initial 
data 



«|*=0 =15 
supplemented with the boundary conditions 



d z u h \ z= i = 0, u 3 | 2=1 = 0, 
d z u h \ z=0 = 0, u 3 | 2=0 = 0. 
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Then v = (vh,0), where Uh is a two-dimensional divergence free vector field 
given by 

v h{t) = ^-^2 I ex P[—T^- u h\kh\ 2 t + ix h -k h Jv i l(k,^ y )d^y 
where kh = (k,£ y ) and 
Vh(k,£y) = 7r[ exp{-i(xk + yQ)v° h (x,y)dxdy, V^^eZxR. 

Proof. Let us first prove that the property d z u = is propagated by equation 
(jl.7p . Using the same arguments as Chemin, Desjardins, Gallagher and 
Grenier in [3] for classical rotating fluids, one can introduce some kind of 
Fourier variable with respect to z, denoted by k^. Since equation (jl.7p is 
linear, it can be easily checked that there is no resonance between Fourier 
modes in k%; in other words, since the only Fourier mode at time t = is 
&3 = 0, there is no Fourier mode corresponding to k% ^ for t > 0, which 
means exactly that d z v = 0. 

We infer that for all t > 0, v{t) is a two-dimensional vector field which 
satisfies 

(A -. s d t v h + -L R v h - Vh&hV h = 0, div h v h = 0, 

(4.1) e 

3 Z P = 0, v 3 = 0. 

This leads to 

v h (t) = exp (t (--j 1 + Vh&h^j v h \ t=0 . 

Let us now investigate the precise expression of the operator Lr. First, since 
Vh is divergence free, we have, for all y £ Ft, 



d y / v 2 (-,y) = - / d x vi(-,y) = 0. 
Jt Jt 

Consequently, since Vh 6 L 2 (T x R), 

[ v 2 (t,-,y) = Vt>0, yen. 
Jt 

Taking the x-average of the first component of (|4.ip . we obtain 

dt vi- v h d 2 y I vi = 0. 

JT JT 

This corresponds to the "stationary part" of v e in Theorem 11.31 

Hence Lemma 14.21 is proved for the Fourier modes such that k = 0, where 
k is the Fourier variable associated with x. Thus we now focus on the modes 
such that k ^ 0, or, in other words, on initial data such that J* T v^ = 0. For 
such vector fields, we have, since Vh £ i 2 (T x R) is divergence free, 
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where 

C(t) := rot h v h = d x v 2 - dyV\. 

On the other hand, 

rot h (bv h )) = div h (bv h ) =v h -Vb = (3v 2 . 
Gathering the last two inequalities, we infer that 

d t ( + -d x A~\ - u h A h C = 0. 

e " 

In Fourier space, this leads to 

and thus, setting kh = (k,£ y ), 

ik 

\kh\ 2 \~t\kh 

Since v is a two-dimensional divergence free vector field, for all kh G Z x R, 
we have kh ■ v^(kh) = 0, and thus 

v h {k h )- ^ k h . 

Eventually, we retrieve 



Vh(t,k,£ y ) = --^expli——^t-u h \kh\ 2 t)C\t=o(k,Cy) 



Vh(t, k, £„) = exp (^ | fc | 2 ~ ^ | 2 * - v h \k h \ 2 tj v° h (k, £ y ) Vt, k, £ y . 
Using the Fourier inversion formula, the proof of the Lemma is complete. □ 

5. Three-dimensional propagation 

5.1. Remarks about the qualitative behaviour of three-dimensional 
waves. 

We are now interested in waves having vertical oscillations, that is in the 
solutions to 

d t u + -fiyu L + ( Ylf ) - VhA h u - ed zz u = 0, 
e \71OzPJ 

v ' V • u = 0, 

d z u h \ z= Q = d z u h \ z= \ = 0, u%\ z= q = n 3 | z=1 = 0, 

having zero average with respect to z. 

Once again, we introduce a kind of Fourier variable with respect to z (see 
[1]), denoted by k$, which here is different from zero. The Fourier variable 
associated with the first coordinate x is still denoted by k. 
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If Uh is sufficiently small, we then expect the main dynamics to be given 
by the Poincare propagation operator 

(5.2) L P u = 0yu^ + (f d h ^ 

where p is such that both the incompressibility constraint and the boundary 
condition are satisfied. 

• A very rough analysis shows that fast oscillations with respect to 
y should appear for times greater than e. Indeed, as long as the solution 

(u,p) to 

edtu + Lpu = 
depends slowly on y, the pressure which satisfies 

III 

-{d X x + dyyp + -^d zz )p = — /3yd x u 2 + -d y {(5yui) 

can be approximated in the following way 

V = A {-ikf3yu 2 + d y {(3yui)) = O(e). 

In particular, at leading order, the singular penalization behaves as in the 
compressible case 

(L P u) h ~ (3yuj^ 
Plugging this Ansatz in the evolution equation leads to 




which is relevant only for very small times, but indicates that a fast depen- 
dence with respect to y can be expected. 

• On the other hand, we do not expect (u,p) to behave as a function 
of y/e only. Such a property, together with usual integrability conditions, 
would indeed imply that the solution (u,p) concentrates on small times in 
the vicinity of y = 0. As previously, a rough analysis based on the change 
of variable Y = y/e and on some asymptotic expansion of Lpu 

Lp~u ~ (o,ik 3 (k 2 3 - dYYrHPYdyui)) 

shows that "concentrated functions" are not stable under the penalization 
hp. 

The mechanism we want to study involves therefore both scales y 

and y/e, and results from a balance between rotation and vertical oscilla- 
tions, which is the main novelty here. Note indeed that previous works on 
rotating fluids consider either the case when the effect of rotation is domi- 
nating (macroscopic layer of fluid) [4] or the case when vertical oscillations 
hold on very small scales and can be averaged (shallow water approximation) 

Semiclassical analysis seems therefore to be the relevant tool to study this 
problem, insofar as it allows to separate both scales in a systematic way. 



/3-PLANE MODEL FOR ROTATING FLUIDS IN A THIN LAYER 



35 



• Note finally that, if the horizontal viscosity is such that 3> e 2 , then 
because of the small scale in y, we expect all the energy to be dissipated 
on a small time interval, leading to some boundary layer effect (see the 
discussion in paragraph 15. 5ft . 

In order to exhibit a non trivial propagation, we will assume in all the 
sequel that 

Vh = o(e 2 ). 

We therefore start with the study of the 3D propagation without dissipation. 
We will then check a posteriori that the viscous dissipation introduces only 
small error terms for any finite time. 

5.2. Semiclassical analysis of the three-dimensional propagation. 

In order to study the propagation of energy by 3D waves, a natural idea is 
then to get a polarization of Poincare waves, i.e. to obtain a diagonalization 
of the system 

edfU + Lpu = 

in the limit e — » 0. We first use the incompressibility constraint to rewrite 
the propagator in the form of a 2 x 2 matrix of pseudo-differential operators. 
We indeed have 

-A e p := -e 2 (A h + -^d zz )p = -e(5yd x u 2 + ed y {(5yui) 

from which we deduce that 

eB u, 4- ( -^d y A 7 H(3y.) -dy ■ -eW xx A 7 HM\ Uh _ Q 
ed t u h + ^ y . _ € 2 d 2 yA7 i m e 2 dxdyA7 i iM )u h -0. 

Our first goal is then to perform a suitable change of variables leading to 

£dtV+ (HHd x ,dz,y,ed y ) ) v = 0{n . 

\ H e (d x ,d z ,y,edy)J 

In all the sequel, for the sake of simplicity, we will consider a single Fourier 
mode in (x,z), and denote by (k, £3) G Z x Z* the associated wavenumber. 
Any solution is indeed a superposition of such waves. We will denote abu- 
sively H^(k, ks,y, ed y ) the Fourier transform of H 7 1 {d x , d z ,y, ed y ). 

We are then brought back to study the propagation of waves by the 
scalar pseudo-differential operator H^(k,k3,y,ed y ), which can be done for 
instance using classical results on the Wigner transform. For such scalar 
skew-symmetric pseudo-differential operators, we indeed know [12] that en- 
ergy is propagated according to the hamiltonian transport equations 

d t f + {h ± j} = o, 

where h ± (k, k%, y, £) is the semiclassical principal symbol of Hf^{k, k%,y, ed y ). 

Note that the time scale over which one has a macroscopic propagation 
of the energy is inversely proportional to the size of the oscillations. Such a 
property can be seen very simply on equations with constant coefficients 

edfV + h(ed y )v = 
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remarking that the group velocity 

has a finite limit as ek 2 — ► £. 

What we are finally able to establish is the following Proposition 

Proposition 5.1. Let vP £ L 2 (lu) be a compactly supported divergence free 
vector field such that J vPdz = , and let u £ C(R+,L 2 (o;)) 6e i/ie solution 
of equation (|5.ip mtt initial data u° . 

Then the L 2 norm of Uh(t) on any fixed compact converges to as t — > oo. 

In other words, 3D waves are dispersive, but only on times of order 1. 
Note that, in the case of a macroscopic layer of fluid, the velocity group 
of Poincare waves is much larger (typically of order 1/e); see for instance 

Furthermore the vertical component U3 of the velocity will not remain 
bounded, as is usually claimed in formal derivations leading to shallow water 
models. 

5.3. Reduction to a scalar situation. 

The first step of the proof follows a method initiated in [5]. 
• We first compute a kind of characteristic polynomial for the matrix 
of pseudo-differential operators 

e(-A e ) 2 )-Hked y (py) -py + e{- A e )~ l k 2 f3y> 

A simple way to obtain a scalar equation is to proceed by linear combination 
and substitution. 

Because the solution is expected to depend both on y and y/e (whatever 
the initial data), ed y is a O(l) operator like multiplication by any function 
of y. We then apply usual rules of semiclassical analysis : 

6^ = 0(1), y = 0(l), 

and any commutator has smaller order 

[ed y ,y] = 0(e). 

Keeping only leading order terms, we get 

irux - (3yu 2 = 0(e), 
Pyui + ed y (k 2 - (ed y ) 2 )~ l ed y (l3yui) + iru 2 = 0(e) 

so that 

PVm + ed y (k 2 - (edyfy'edyiPV^) - t 2 u 2 = O(e) 
or equivalently 

(5.3) k 2 ((3y) 2 u 2 - r 2 (k 2 - (ed y ) 2 )u 2 = 0(e) 
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since commutators provide higher order terms with respect to e. Note that 
one can also compute an exact pseudodifferential relation (which is actually 
a polynomial of degree 6 with respect to r) by keeping all the terms 

(5.4) P(e,y,ed y ,r)u 2 = 0. 

Note that, contrarily to [5], as we will only consider times of order 1, we 
do not need to compute subsymbols, so that we could also proceed directly 
using symbolic calculation and diagonalize the matrix 

/ -0y\ 

Anyway, we expect the roots to the following polynomial to play a special 
role in the propagation : 

(5.5) P(0, y, £, r) = k 2 ((3y) 2 - (fc§ + f)r 2 



• We can actually prove that there exist pseudo-differential operators Hf 1 
with principal symbols 



h 



± 



such that edt/J.^ = iH^/j.^ implies that 

v ± := Id W t ± satisfies flED up to 0(e°°), 

where are scalar functions. 

This result is actually a variant of the main Lemma in [5j. (Indeed the 
exact dispersion relation depends here explicitly on e.) 

Lemma 5.2. |5j Let P t = P(e, y, £, r) be a smooth function such that 
9 r P |P=o 0? anc ^ let h = h(y,£) be any continuous root of 

P(0,y,C,h(y,0) =0. 

Then there exists a pseudo-differential operator H e = H e (y, —ied y ) with 
principal symbol h(y, £) such that: 

(5.6) H e ip = TTp =^ P e>T ^ = 0(e°°) 

where P €jr is a pseudo-differential operator of full symbol P(e,y,^,r). 

The proof of this lemma relies on pseudo-differential functional calculus, 
and uses various quantifications to make the computations as simple as 
possible. For the sake of completeness, we recall here the main arguments, 
but refer to [5] for details. 
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At first order, we have 



p 



e 

I J„J A„JI 



So the principal symbol of P £)T is P(0, y, £, h(y,£)) which, by assumption, is 
0. 

For the e°° result, it is enough to repeat the same argument with h e ~ 
h + Yl ek hk- We obtain 

P(e, y, £, fr e ) + ^ e*Q fc (/i, . . . , ^<9f fc £ ) = 0, 
k>i 

that can be solved recursively under the condition <9 T P |p =0 / 0. 

• We further obtain a decomposition of any initial data on the eigenstates 
of the scalar propagators H^ 1 . 

For all u°, there exist such that: 

UQh = (~ ( Py " edy^edyiPv))' 1 W - ed y k 7 \iekPy)\ 

j \ M / 

+0(e°°) 

=: ^Q^ + 0(e 00 ). 
± 

where A e := e 2 <9 2 y — e 2 /c 2 — k\. The vertical component is then entirely 
determined by the divergence-free condition. 

To prove this result, one first remarks that the leading order symbol of 
the matrix namely 

\k 3 \sgn( y ) IMsgnfe) 



1 1 



is invertible. 



The inversion of the matrix (Q+ Q e ) can then be done symbolically at 
any order. 



5.4. Dispersion of energy. 

Standard arguments of semiclassical analysis allow then to control the 
propagation of energy for the scalar equations 

edtnf + iHfuf = 
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• Because Lp is skew-symmetric (in some weighted L 2 -space), we have a 
uniform control on the L 2 norm of Uh 

IKIl£» + e 2 |H|£2 (u;) = K|| 2 2(£j) + e 2 lh3lli» 

These uniform a priori estimates allow to establish the convergence of the 
remainders in the equations for the Wigner transforms 

f?(t,y,Q := ^1 e^y'ntty-ey'Y^y + ey'W 
We therefore have 

d t ft + {h ± jt) = 0{e). 
For detailed computations leading to that estimate, we refer for instance 
to Cm or [12 : 



Lemma 5.3. [T2] Let // ± be any fixed function of L 2 (non-oscillatory) . 
Assume that 

• iHf is self-adjoint on L 2 ; 

• there exists a G R such that is of order a uniformly as e — > 0; 

• the Weyl symbol of satisfies 

hf = h + ehf + o(e) uniformly in Cf£ c . 

Then the Wigner transform /<r(i, y, £) of [if{t) converges locally uni- 
formly in t to the continuously t-dependent positive mesure f^, solution 
to 

d t f ± + {h ± ,f ± } = 0. 

In other words, the energy associated to the ± mode is transported along 
the characteristics of the hamiltonian /i ± : 

(5 7) dt * 

dt dy 1 ^ >' 

• The previous ID hamiltonian systems are of course integrable. The 
bicharacteristics are indeed included in the level lines of h , which are hy- 
perbola as shown in Figure 2. 

A rapid inspection of the large time asymptotics show that trajectories 
cannot be trapped in some compact. This would indeed imply that there 
exists either some stationary point or some turning point. But H(i) is a 
monotonic function 

dE^ (3k 3 sgn(Y(t)) 
dt T ytk'i + Z 2 ^) 
which converges necessarily to infinity. 

For any fixed compact, we can even get an explicit estimate of the exit 
time since 

|3(t)-3o| >(3t, 
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Figure 2. Bicharacterictics associated to Poincare waves 



/3-PLANE MODEL FOR ROTATING FLUIDS IN A THIN LAYER 



11 



from which one deduces a similar estimate for Y(t) : 



Y(t)\ > -^^l + (/3t-|H |) 2 



Note that, since the initial data uq we consider is supposed to depend only on 
the slow variable y, all bicharacteristics we are interested in satisfy Ho = 0. 
By definition of the wavefront set, we finally obtain Proposition 15.11 

Remark 5.4. The qualitative behaviour of Rossby and Poincare waves ob- 
tained here, i.e. in the case of a thin layer of fluid with rigid lid, is very 
different from the one exhibited in shallow water approximations (see [5] ). 
Note that, in both cases, Rossby waves are easily identified because they are 
directly linked to the inhomogeneity of the Coriolis force, in particular they 
always propagate eastwards. 

Here the energy associated to Poincare waves propagates much slower than 
the energy associated to Rossby waves. The point is that fast oscillations with 
respect to latitude y, which are generated spontaneously for vertical modes 
but not for purely 2D Rossby waves, slow down the propagation. Maybe it 
would be physically relevant to consider initial data that depend already on 
the fast variable y/e. 

The other point which should be discussed is the influence of the free- 
surface. But, at the present time, we have no convenient mathematical tool 
to study the propagation of waves in such a complex geometry. 

5.5. Influence of the viscosity. In the case when = o(e 2 ), an easy 
computation based on the energy estimate shows that the viscous dissipation 
does not modify the propagation for finite times. 

More generally, we could extend the previous study considering the whole 
viscous Poincare propagation operator 



where p is such that both the incompressibility constraint and the boundary 
condition are satisfied. 

The diagonalization process is of course unchanged since the dissipation 
operator is scalar. The only difference is therefore that one has now to 
control the propagation of energy for the scalar equations 



A standard computation (reported for instance in Proposition 1.8 of |12j ) 
shows that the Wigner transform then satisfies the following damped trans- 
port equation 



(Note that the symmetric part of the propagator occurs at leading order in 
e, which can be seen by easy symmetry considerations.) 





edtut + iHfuf - Vh^h^t ~ td zz fif = 0. 



dtft + ^\i\ 2 ft + {h ± jt) = o{l). 



12 
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We then deduce that 

• if Vh *C e 2 , the energy is propagated according to the bicharacteris- 
tics associated to h^, as stated in Proposition 15.11 ; 

• if Uh 3> e 2 , the energy contained initially in the Poincare modes is 
dissipated on a very short time, leading to some initial layer phe- 
nomenon ; 

• if Vh ~ e 2 , the dynamics is a combination of both phenomena, as 
shown by Duhamel's formula 

Y ± (t, y, 0,^((t, y, 0) exp (4^|H(t, y, £)| 2 t) = f°(t, y, . 

Note in particular that the energy associated to Poincare modes has a 
super exponential decay, since |E(t, y, £)| — ► oo along any trajectory. 



6. Derivation of the thermocline 

This section is devoted to the proof of Proposition 11.61 which relies on 
classical elliptic arguments. The main difficulty lies in the fact that the 
equation on 9 is degenerate in the horizontal variables. We first prove the 
existence of 9, along with some H 1 estimates, and then we prove the con- 
vergence. 

Throughout the proof, we assume that the wind stress a vanishes at 
sufficiently high order near y = 0, so that there is no need for a truncation 
(see section [2]) and the function u stat does not have any singularity. 

• A priori estimates on the function 9: 

Let 9 £ L 2 (w/ l ,iJ 1 ([0, 1])) be any solution of (fTTTjh Multiplying (11171) 
by 9 and integrating on u, we obtain 

A / \d z 9\ 2 = -~ / u int -nj 2 + \l ^rv. i'V. i- A/ 0-A-. A 



(6.1) -\ [ v%f z=l 9 2 1 + \ [ 9 1 d z 9 {z=1 . 

According to section [3l we have 



Ant _ ^2 B Ol 

M 3| 2 =i- b Oyy 



We assume that a is such that the right-hand side belongs to L 00 (cj/ 1 ). We 
now evaluate d z 9\ z= i : we have 



l pi 



Xd z 9 {z=1 = XI d zz 9= I u int ■ V9 



o 



(6.2) = div,«* / 9) + uf z= A- 
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Recall that u^, defined in section El is independent of z, while is linear 
with respect to z. Consequently, the function 9 depends on and z, and 

xf 9 x d z 9\ z=x = [ 9i(^v h (yt f S) + t^ =1 1 ) 

J ujh J Uh V J / 

(6.3) = - f 9ul nt • V fc fli + / vftUfi. 

Using the identity 

0(.,z) =0 X - I dj(-,z') dz', 

J z 

we deduce that 

(6-4) l|0||£ 2 (u/) < II^iIIl 2 ^) + II^^IIl 2 ^)- 

Gathering (|6.ip . (|6.3p and (|6.4p . we infer that 

a/|M 2 = 5/ 4?= A 2 - / • v ft 0i 

<r ll„,* n * II 1 1 /a ||2 

+ K"'IM|V/AIMII#ilb + \\d z 9\\L*)- 
Using the Cauchy-Schwarz inequality, we obtain eventually 

(6.5) a ! \d z e\ 2 < IKtilk-ll^ill^ + lk^lli-llv^ill^ll^lU^ 

J UJ 

1 1 ||„,int||2 11x7 a i|2 

Inequalities (|6.5p and (|6.4p entail that any solution 6 of (|1.17|) is bounded 
in L 2 (iOh, if 1 ([0, 1])) by a constant depending only on A, 9\ and u int . 

We now derive estimates on the horizontal derivatives in a similar fashion: 
we have 

(6.6) - \d zz v h e + (u mt • v) v h e = - ( v fc <* ) • v h - v h v$*d t S. 

Multiplying the above equation by V/j# and integrating by parts, we have, 
using the boundary conditions, 

- / a**v fc £ • V h 9 = [ \d z V h 9\ 2 - I d z v h 9 lz=1 ■ V h 9i 

J UJ J UJ J UJfr 

d z V h 9\ 2 + ! aJ, z=1 A^i. 

Jui h 

Using equation (|6.2p . we express d z 9\ z= i in terms of 9 and #i. Integrating 
by parts once again leads to 



d z 9 {z=1 A h 9x 



< \ {hh%~\\9h40i\\ H 3 + ||4 n 'lk°°l|0ill^lNM • 
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On the other hand, since u l ^ =0 = 0, we have 



2 ! v)v h o\-v h e = [ uf z=1 \v h e lz=1 \ 2 - [ u^ z=0 \v h 9 lz=0 \ 2 

= [ 4tilv,Al 2 . 

Ju) h 

The two terms in the right-hand side of (|6.6p can easily be evaluated in L 2 
using the estimate on d z 9; there remains 

a / \d g v h e\ 2 <c+\\v h u%*\\ L °°\\Vh^h> 

where the constant C depends on A, ||it mt ||i,°° and ||^i||jy3. 
Assume that 

this assumption is discussed in Remark 11.71 following Proposition 11,61 Then 



/ |3 2 v^| 2 + / \v h e\ 2 < c, 



where the constant C depends on A, 9\ and u int . These estimates easily lead 
to the existence of a solution 9 of equation (11. 17ft ; the uniqueness of 9 follows 
from the estimates above with 9\ = 0. The same method also shows that 
under condition (|1.14D on V^"*, Dffi is bounded in L 2 (uj h , H 1 ^, 1]). Plug- 
ging this estimate back into (16.6|) . we deduce that V^(9 G L 2 (ujh, H 2 [0, 1]), 
and thus that V h 9 is bounded in L 2 (uj h , W 1 ' 00 ^, 1])). 

Concerning the function 9 BL , the existence and uniqueness are obvious; 
we have merely 

(xk.O-^V^-^iw ) (A±(xfc))3 . 
• Proof of convergence: 

We construct an approximate solution of (|1.15p as follows: we set 



x h , z) = 9{x h ,z) + e9 BL ( x/j, - — - ) + e0(s ft , z 



1 - z 



where the function 9 is defined by 



1 

z-1) 

Notice that by construction 



9(x h , z) = (z- l)-^d c 9^(x h ) - 9f ( i (x h ). 



d z 9? pp n = 0, 9* p \ 

* \z=0 ' z=l 



'1- 



Moreover, using the definition on 9 BL , it is easily proved that 9 = O(l) in 
W 2 >°°(u) (provided the stress a is smooth and vanishes at a sufficiently high 
order near y = 0). 
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Consequently, 

(6.7) -\d zz e app - \e 2 A h e app + u stat ■ ve app 

= u BL ■ V h {9 - 61) - Xe 2 A h 9 + u§ L d z 9 + v int ■ V9 
-\e 3 A h 9 BL ( x h , —) + eu stat • V9 BL ( x h , — 



€ 

-Xe 3 A h 9 + eu stat ■ V0. 

According to the results of sections [2] and [3l we have 

11^11^= O(^), \\v int \\ L i = o(e), 

4< tat h™,\W3 at h°° = 0(1), e\\u stat \\ L 2 = 0(V~e). 

These estimates, together with the ones derived above on 6, enable us to 
bound all the terms in the right-hand side of (|6,7p . except for the first one. 
Using Hardy's inequality, we have 



W BL -V h ( 



1l )\\L 2 (w) 



< \\( l - z ) u h L \\ L ~ iuih , L 2 m])) \\0--*) 1v ^(^-^)|Il2 ( ^,l-([o,i])) 

< cv~e Wd^nie-e,) || L2Kjioo([0>1])) . 

Thus 6 app is an approximate solution of (11.150 . with an error term o(l) in 
L 2 (u>). As a consequence, 9 — 9 app satisfies 

-\d zz {9 - 9 app ) - \e 2 A h {9 - 9 app ) + u stat ■ V(0 - 9 app ) = o(l), 
(9-9 ap % =1 = 0, d z (9-9 app ) lz=0 = 0. 

Multiplying the above equation by 9—9 app and using the Poincare inequality, 
we prove that 

\\d z (9-9 app )\\ L2 ^=o(l\ 
and thus the Proposition is proved. 
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